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ABSTRACT 

We show that the luminosity of a star forming galaxy is capped by the production and subsequent 
expulsion of cosmic rays from its interstellar medium. By defining an Eddington luminosity in cosmic 
rays, we show that the star formation rate of a given galaxy is limited by its mass content and the 
cosmic ray mean free path. When the cosmic ray luminosity and pressure reaches a critical value as a 
result of vigorous star formation, hydrostatic balance is lost, a galactic scale cosmic ray-driven wind 
develops, and star formation is choked off. Cosmic ray pressure-driven winds are likely to produce 
wind velocities in proportion to and significantly in excess of the galactic escape velocity. 

It is possible that cosmic ray feedback results in the Faber- Jackson relation for a plausible set of input 
parameters that describe cosmic ray production and transport, which are calibrated by observations 
of the Milky Way's interstellar cosmic rays as well as other galaxies. 

Subject headings: galaxies: formation - galaxies: fundamental parameters - galaxies: starburst 



1. INTRODUCTION 

In the context of galaxy formation, "feedback" pro- 
cesses are commonly invoked to limit the formation of 
stars. Feedback can be "intrinsic" or "extrinsic," where 
an intrinsic feedback mechanism regulates star forma- 
tion directly at the molecular cloud level and an extrinsic 
mechanism controls the supply of stellar fuel on galactic 
scales. 

There is long-standing observational evidence that ex- 
trinsic feedback mechanisms operate. For example, the 
luminosity of local galaxies correlates with the stellar ve- 
locity dispersion, and thus gravitational binding energy, 
of the entire galaxy itself (Faber & Jackson 1976). Two 
popular extrinsic feedback mechanism arc thermally- 
driven winds powered by core-collapse SNe (see e.g., 
Chevalier & Clegg 1985; Leitherer et al. 1999) and 
momentum-driven winds powered by starburst radiation 
(Scoville 2003, hereafter SCO; Murray et al. 2005, here- 
after MQT). In this work, we consider another possibil- 
ity: galactic winds driven by cosmic rays. 

1.1. Basic Idea: The Eddington Limit in Cosmic Rays 

There is a theoretical upper-limit for the luminosity of 
gravitationally bound systems, the Eddington limit L c( jd; 
given by 

ATrGM cnc m p c 

-^edd = U] 

0~T 

where M enc is the enclosed mass and <7t is the Thom- 
son cross section. Quasars, the most steadily luminous 
objects in the Universe, are thought to be powered by 
accretion onto a black hole and thus their maximum lu- 
minosity is given by 

lO^Mgergs" 1 (2) 
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where Mg is the quasar black hole mass in units of 
10 9 Mq. In the case of elliptical galaxies, the mass of 
the black hole's host galaxy is ~ 10 2 — 10 3 times that of 
the black hole (Magorrian et al. 1998). However, rather 
than being the brightest persistent sources in the Uni- 
verse, typical luminosities of bright star forming galaxies 
range from 10 45 — 10 46 ergs -1 - at least a factor of ten 
less powerful than their quasar counterparts. 

Photons are not the only type of radiation produced 
from the act of star formation. For every core-collapse 
supernova, a fraction e ~ 0.1 of the released energy 
~ 10 50 erg is converted into cosmic rays (CRs). Unlike 
photons, cosmic rays interact many times with the mag- 
netized interstellar medium, and exchange momentum 
with it, through pitch- angle scatterings with Alfvcnic dis- 
turbances. Though the energy injection rate of CRs is 
insignificant in comparison to other by-products of star 
formation such as photons, winds, and supernova shocks, 
the rate at which they inject momentum is not. 

In the case of the Milky Way, the CR diffusion time 
idiff.CR ~ 10 7 yrs over a characteristic scale of H ~ lkpc. 
Thus the characteristic mean free path for a Galactic CR 
is given by 
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for 1-10 GeV CRs. Whereas the Thomson mean free path 
At ~ 10 24 cm. Thus, CRs are much more highly coupled 
to the Galactic matter in comparison to photons. 

Due to the higher level coupling, the threshold for the 
maximum allowable CR flux and luminosity is lower in 
comparison to the electron scattering case. A rough theo- 
retical upper limit for the CR luminosity £ c dd cr is given 
by 
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Ultimately, the ratio A CR /At may be thought of as the 
inverse of the ratio of the effective cross section for a CR 
to interact with a unit of plasma mass divided by the 
corresponding value for a photon. 
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1.2. Plan of this work 

In the next section we present the main results of this 
work. First, a brief overview of the observed properties of 
bright starbursting galaxies is given in H2.ll The basics 
of a momentum-driven CR feedback mechanism is given 
in £12.21 In £12.31 and to some degree in <j2]4] as well, we 
weigh the advantages and disadvantages of the various 
galactic feedback mechanisms and then compare them 
to our CR feedback mechanism. In $3] we discuss some 
of the relevant microphysical properties of CRs that are 
relevant for CR feedback and in §3.5.1l we speculate as to 
whether or not CR feedback is responsible for the Faber- 
Jackson (1976) relation. We conclude in |j4j 

2. COSMIC RAY FEEDBACK IN STARBURSTING 
GALAXIES 

Galactic scale CR-drivcn winds have been studied from 
a theoretical point of view for many years (see e.g. 
Ipavich 1975; Breitschwerdt et al. 1991). However to 
our knowledge, the connection between star formation 
and wind power as a function environment was not con- 
sidered. 

Recently, CR feedback in galaxy formation has been 
considered by Jubelgas et al. (2006). They found that 
CR feedback can quench star formation for low mass 
galaxies with velocity dispersions a < 80 km/s. In what 
follows, we argue that CR feedback potentially chokes 
off star formation even for galaxies that are much larger, 
where the velocity dispersion a can reach values up to 
~ 300 km/s. We believe that our conclusion differs from 
theirs due to a difference in the handling of CR transport, 
a point that is discussed in H3.4I 

Before we delve into the physics of CR feedback and 
transport, some of the observed properties of starburst- 
ing galaxies are summarized. 

2.1. Observed Properties of Starbursting Galaxies 

As previously mentioned, galaxies and even starburst- 
ing galaxies are dim in comparison to their quasar coun- 
terparts. Bright young and gas rich starbursters possess 
characteristic bolometric luminosities of Lboi ~ 10 45 — 
10 46 ergs" 1 at the massive end. Starbursters are powered 
by UV photons emanating from the surface of massive 
stars in which the overwhelming majority (~ 90 — 99%) 
of the energy release is reprocessed into the infrared, a 
result of UV opacity onto dust grains. 

Although starbursting galaxies exhibit rather complex 
structure and characteristic properties may vary greatly 
with location, some representative values of the densi- 
ties, outflow velocities, and length scales can be obtained 
from observations. Typical H2 number densities for ex- 
treme starburst regions vary between 10 2 - 10 4 cm . 
These densities are averages over regions which extend 
~ 100—500 pc from the nucleus (e.g., Downes & Solomon 
1998). High velocity winds are a generic feature of star- 
bursts, both locally and at high redshift (Veilleux et al., 
2005). Outflow speeds in excess of the escape velocity 
are observed irrespective of the mass of the host (see e.g., 
Fig. 7 of Martin 2005). Mass outflow rates can also be 
estimated, although less robustly. Outflow rates are typ- 
ically inferred to be anywhere from ~ a few to ~ 100% 
of the star formation rate (Veilleux et al. 2005; Martin 
2005). These estimates are sensitive to the length scale 



used to convert the observed gas column to the mass 
outflow rate, which is difficult to determine. 

There is also copious evidence for high velocity winds 
in the spectra of several z > 1 galaxies. This evidence 
comes from blue-shifted absorption lines and P Cygni 
Lya emission line profiles in the rest frame UV spectra 
(Veilleux et al. 2005). Among these high redshift galax- 
ies, the lensed source, MS 1512-cB58, has the highest 
quality spectra. The bulk outflow velocity is 255 km s _1 , 
with absorption extending to ~ 750 km s _1 , well above 
the escape velocity (Pettini et al. 2002). These outflow 
velocities are not uncommon, and are, in fact, typical of 
luminous Lyman break galaxies (Adclbcrgcr et al. 2003, 
Shapley et al. 2003). 

2.2. Maximum Luminosity of Starbursting Galaxies 

As previously mentioned, CRs arc an energetically sub- 
dominant by-product of star formation. But, due to their 
small mean free path with the magnetized interstellar 
medium, the collective momentum response with a given 
galaxy's gaseous component might have dramatic conse- 
quences. In order to determine under what conditions 
CR production will heavily influence the dynamics of a 
galaxy, we consider galactic hydrostatic balance were the 
CR force is balanced by gravity alone. 

For the sake of simplicity and a lack of theoretical and 
observational constraints, we take a one-zone approach. 
That is, we collapse all of the information regarding the 
structure of a given galaxy into a single characteristic 
value taken at a characteristic length scale. Later on, in 
£12.3l we justify this approach. Furthermore, from here on, 
we assume that the mass distribution of a given galaxy 
resembles that of a singular isothermal sphere whose rel- 
evant properties are briefly listed in Appendix [X] 

In the interest of clarity, we derive the maximum lu- 
minosity in CRs - and incidentally photons as well - 
of young star forming galaxies by examining hydrostatic 
balance in three different ways. 

2.2.1. Hydrostatic Balance 

In order for the CRs to limit the inward flow of gas, 
thus quenching star formation, they are most likely re- 
sponsible for the bulk of the galactic pressure. In the 
diffusion limit, the CR pressure is related to the CR flux 
by the relation 

CR — CR ! V u 7 

where r CR is the CR proton optical depth. The CR en- 
ergy flux F CR is related to the star formation rate per 
unit area (i SF by the relation 

F cr - CcrAsfC 2 - (6) 

In the above expression e CR and (i SF are the efficiency 
of converting mass into CR proton energy and the star- 
formation rate per unit area, respectively. If we assume 
the galaxy in question resembles an isothermal sphere, 
the flux of CRs can be related to a star formation rate 
within a given radius 

F CR ~ 6 x lO-^gm^^crgcm^s- 1 (7) 

where eg, m SF , -Rkpc is the efficiency of converting rest 
mass into CR energy in units of 10 -6 , star formation rate 
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in units of solar masses per year, and galactic radius in 
units of kpc, respectively. 
Now, we may write the CR pressure as 

P CR ~ 3 x 10- 13 r CR e 6 m SF i?- p 2 c ergcm- 3 . (8) 

In our feedback picture, the outward diffusion of CRs 
halts gas accretion and thus star formation once CR pres- 
sure is almost entirely supporting the gaseous component 
of the galaxy against the inward force of gravity. In that 
case, hydrostatic balance dictates 

Pen ^P^^X 10-7^4)0 (9) 

for an isothermal sphere. With this, a criteria relating 
star formation rate and CR optical depth to the mass 
of the galaxy, or equivalently, its velocity dispersion is 
given by 

m SF ~10 3 T 3 - 1 e6 -V 9 o. 1 ^oo. (10) 
which docs not explicitly depend upon galactic radius. 
Note that T3 and f go 1 are the CR optical depth in units 
of 10 3 and the gas fraction normalized to 0.1, respec- 
tively. The value for m SF in the above expression seems 
enormous. Interestingly, this is truly an upper limit. 

First, it is unlikely that CR feedback can affect a galaxy 
whose velocity dispersion is much larger than 0300 since 
only ~ IkeV per baryon (which corresponds to a velocity 
dispersion a ~ 300 kms -1 ) is available in the form of CRs 
for an efficiency eq. From the arguments in Appendix [Dl 
it follows that the maximum value of the optical depth 
for a wind with a terminal velocity v w ~ 0300 is of order 
r CR ~ 10 3 . Furthermore, a galaxy with a velocity disper- 
sion a ~ 300 kms certainly represents the high end of 
galactic mass and gravitational potential and thus, one 
would not expect values of the star formation rate in 
excess of eq. (fT0| . 

The result given by eq. (|10p yields a maximum star 
formation rate that is nearly identical to the value ob- 
tained by eq. (19) of MQT as long as r CR and e CR are 
independent of galactic mass and velocity dispersion. In 
what follows, we explore the meaning of this apparent 
coincidence. 

2.2.2. (Cosmic) Radiation Pressure-Driven Wind 

Instead of balancing CR pressure with gravity, we now 
phrase hydrostatic balance in terms of the CR luminosity 
and enclosed mass. Again, by assuming spherical sym- 
metry, hydrostatic balance is given by 

«CR„ KCR L CR GM(r) 

^cr = — 2 ~ 2 — ■ ( LL ) 

c c 47rr z r A 

If we integrate over mass i.e., dm = pd 3 x = 4irpr 2 dr, 
then opacity may be replaced with optical depth. We 
have (see e.g., Lamers & Cassinclli 1999) 

r CR _ f g G M 2 (R) 

L ™—- R? (12) 

If the mass distribution of the galaxy in question resem- 
bles an isothermal sphere, then we may write 

L CR r CR * (13) 

which is strikingly similar to the Faber- Jackson (1976) 
relation. In terms of typical scalings for galaxies 

ic R ^cR - 3 x 10 46 jgo.^oocrgs- 1 (14) 



where f go 1 and 0200 are the gas fraction and velocity 
dispersion normalized to value of 0.1 and 200 kms" 1 , re- 
spectively. Note that eq. (fT4")) is identical to eq. (fTTTj) . 

For every gram of matter that is assembled into stars, 
CRs are responsible for only ~ 10~ 3 of the energy release 
while photons are responsible for the overwhelming ma- 
jority. Consequently, the CR luminosity is smaller than 
the starburst photon luminosity by a factor of ~ 10 3 . 
Motivated by the assumption that every UV photon re- 
leased from the surface of a massive star is absorbed 
by dust grains and subsequently downgraded into non- 
interacting infrared radiation, MQT employ the "single 
scattering approximation," which sets the photon opti- 
cal depth r 7 — 1 for all galaxies. If t cr ~ 10 3 as in the 
case of the Milky Way, then the maximum star formation 
rate, CR luminosity, and photon luminosity given by eq. 
pT|) is equal to the upper limit derived by MQT in the 
case where photon feedback determines the maximum 
galactic luminosity. 

Note, that we have yet to consider the dependencies of 
bulk galactic properties such as the velocity dispersion 
a and gas fraction f g upon the optical depth r CR and 
efficiency e CR . 

2.2.3. The Eddington Limit in Cosmic Rays 

Now, we exchange CR optical depth t cr with CR opac- 
ity kcr in order to obtain an equivalent expression for 
the maximum CR luminosity L cn . In order to do so, we 
must determine the effective cross section or mean free 
path that a typical CR has with magnetized interstellar 
matter. 

Let us take the 1-10 GeV CR mean free path to be 
A CR ~ 1 pc, approximately the Milky Way value, as a 
benchmark. In order to obtain an Eddington limit, we 
compare the CR mean free path A CR to the Thomson 
mean free path At i-e., 

£cdd,CR - L cdd ^- = AttGcM (r)p(r) A CR (15) 
at 

or 

Lcdd.cR * 3 x lO^AJ^.^oo-Rkpc ergs" 1 (16) 

where X 1 is the CR mean free path in units of a parsec. 

The arguments surrounding eq. (|15[) are identi- 
cal to the Eddington argument put forth by Scov- 
ille (2003, hereafter SCO) regarding UV photon opac- 
ity on dust. His primary motivation was to ascertain 
the supposed common origin of the light to mass ratio 
L/M ~ 5OOL0/M0 of both self-gravitating giant molec- 
ular clouds in M51 and the central region of the dust-rich 
luminous starburst galaxy Arp 220. 

We define a CR to mass ratio TZcr = L CR /M(r), which 
for an isothermal sphere is given by 

K C R ^ 2/ g CT200^kp 2 A^O' ( 17 ) 

where TZq — Lq/Mq ~ 2ergs _1 g _1 is the light to mass 
ratio of the Sun. If the efficiency of generating photon 
power from the act of star formation is e S F/ e CR ~ 10 3 ; 
then the upper limit for the light to mass ratio of star- 
bursting environment, TZsf, is given by 

ft SF ~ ^Tlcn ~ 2 x 10 3 f g al m R^ c \n & . (18) 

CR 
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The gas fraction is roughly f g ~ 1/2 for both sets of 
environments considered by SCO. For a CR mean free 
path A CR ~ lpc, the above value for TZsf is reasonably 
close to the measured value for Arp 220, given that the 
velocity dispersion of Arp 220 is slightly different than 
C200 1 the mass distribution is not perfectly described as 
an isothermal sphere, and the observed Lsf (or the bolo- 
metric luminosity) is only accurate to a factor of 2 or 
so. Giant molecular clouds have typical dimensions of 
~ 10 — 30 pc and velocity dispersion a < lOkms -1 . 
Therefore, according to the form of 7£sf given by eq. 
(fT8|) . both systems will be limited to a similar light to 
mass ratio resulting from CR feedback. In Appendix [El 
we discuss the viability of intrinsic CR feedback in giant 
molecular clouds. 

2.3. Cosmic Rays vs. Photons and Supernovae 

Now, we assess how the three major by-products of 
star formation i.e., starlight, supernovae, and cosmic rays 
compare with one another in terms of their ability to in- 
ject momentum into the interstellar medium. Special 
attention is given when comparing CRs and stellar pho- 
tons. Our arguments are given below. 

2.3.1. Global Kinematics of Various Feedback 
Mechanisms 

Energetically, cosmic rays are a highly subdominant 
by-product of star formation. However, if the injection 
of momentum determines whether or not the star form- 
ing constituents irreversibly disrupt the galaxy of their 
birth, then rate of energy injection may not, in fact, be 
determining factor. At the crudest possible level of dis- 
course, the rate at which momentum is deposited in a 
galaxy P is given by 

P ~ E/v ~ L/v (19) 

where v is the characteristic velocity at which the energy 
carrying particles are transported through the interstel- 
lar medium. 

As previously mentioned, L SF > L SN > L CR . However 
at the same time, c > i>sn > ^cr. In the case of corc- 
collapse SNe, 

^sn ~ L SN /w SN ~ 10 33 m SF gems- 2 (20) 

where «sn ~3x 10 s cms -1 is the characteristic velocity 
of a supernova shock. 
For cosmic rays, 

^CR ~ ^Cr/ U CR ~ ^CR T Cr/ C (21) 

~3 x 10 33 e CR tosf gems" 2 , 

which is slightly larger than the value obtained from su- 
pernovae. Here, «cr ~ H/tdis.cR ~ 10 7 cms _1 is the 
characteristic velocity at which cosmic rays diffuse in 
bulk. Note that vcr ~ c/t cr , where t cr is the optical 
depth in cosmic rays given by r CR ~ H/X CB ~3x 10 3 
for the Milky Way. 
In the case of starlight, 

P sf ~ L sr/c~ 5 x 10 32 m SF gems" 2 . (22) 

Therefore all three forms of momentum deposition are 
comparable. If the photons are optically thick, then 
c — > c/t 7 in the above expression where t 7 is the photon 
optical depth. 



2.3.2. Cosmic Rays vs. Photons 

Eqs. (fT0|) . (fT4|) and (fT8|) all arrive at the same conclu- 
sion i.e., for a CR mean path A CR ~ lpc and radiative 
efficiency e CR ~ 10 -6 - values that are reasonably close 
to the Milky Way value - then, the injections and transfer 
of 1-10 GeV CRs will inevitably unbind a given galaxy's 
gaseous component. However, eqs. (|T0|) . (|14|) and (|18|) 
also inform us that for the same choice of parameters, the 
basic observational signature of CR feedback - a maxi- 
mum value for the equivalent physical quantities m SF , 
Lsf and IZsf ~ simultaneously indicate that it is photon- 
feedback, according to the theory of SCO and MQT, that 
caps the luminosity of a given galaxy. 

To our knowledge, there is no simple straightforward 
observational test that distinguishes between the CR 
mechanism presented in this work and the photon-driven 
theory. However, a clear distinction can be drawn on the- 
oretical grounds alone. 

The advantage of utilizing starburst photons as a feed- 
back mechanism, which limits star formation, is that the 
overwhelming majority of the energy and power (~ 99%) 
from the action of assembling gas into stars is released 
in the form of photons. At the same time only ~ 1/1000 
of the energy release and power from star formation 
comes out in the form of CRs. Furthermore, the cou- 
pling between starburst photons and dusty matter is rel- 
atively well understood: UV starlight is absorbed by dust 
grains and immediately downgraded into non-interacting 
infrared radiation. The frequency-weighted opacity for 
this process is roughly a few hundred times the Thom- 
son opacity (SCO; MQT). Finally, measurements of the 
bolomctric photon luminosity of gas-rich starbursting en- 
vironments directly constrain any photon-feedback the- 
ory, whereas for CR feedback a CR proton conversion 
efficiency e CR must be assumed from the starformation 
rate. 

In order to unbind interstellar matter in the largest 
and most massive galaxies, the CR mean free path A CR ~ 
lpc or smaller, which implies that the CR optical depth 
t cr must be large such that r CR ~ 10 3 , or even higher. 
Therefore, any theory of galaxy feedback regulated by 
CRs requires, as in the case of the Milky Way, a large 
CR optical depth. 

We now raise a crucial point: Due to their large op- 
tical depth, each individual CR injects momentum into 
the interstellar medium in the direction opposite a given 
galaxy's gravitational center. Thus, the CR flux and 
pressure gradient do not strongly depend on the spatial 
distribution of injection sites. Furthermore, the fact that 
r CR must be large implies that the one-zone model for- 
mulation invoked in ^2.2.1112.2.31 is well justified at the 
order unity level. 

The momentum recoil of the interstellar medium re- 
sulting from the absorption of UV starlight onto dust 
grains strongly depends on the spatial distribution of 
the photon injection sites. To demonstrate this, we relax 
an important simplifying assumption employed by SCO 
and MQT i.e., approximating that the entirety of a given 
galaxy's UV radiative flux emanates from a single cen- 
trally located point, behind a common dust sublimation 
layer. 

In the Milky Way, for example, there are ~ 10 4 O 
stars (Wood & Churchwell 1989) and for young stellar 
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populations with a Salpeter IMF, O stars dominate the 
radiative power. For the sake of our argument, we as- 
sume that each individual O star is contained within 
a molecular cloud-like environment for the overwhelm- 
ing majority of its life. Note that this is not the case 
in the Milky Way, but is likely to be a good approxi- 
mation in a gas-rich starbursting galaxy. As previously 
mentioned, UV starlight from embedded massive stars is 
reprocessed and downgraded into the infrared by absorp- 
tion onto dust grains (Churchwell 2002). It is this source 
of opacity upon which SCO and MQT build their theory. 
By doing so, they conclude that starburst photons enjoy 
a relatively large momentum coupling with the gaseous 
component of their parent galaxy. 

In Appendix [Bl we briefly outline why the one-zone 
approximation - meaning that the UV photons originate 
from a single point - of SCO and MQT qualitatively and 
quantitatively breaks down on galactic scales. The ob- 
servational consequence is that starlight can only "stir 
up" the gas to characteristic velocities of order v ~ 1/10 
of the velocity dispersion a on length scales of order 
Ax ~ lOpc, respectively. It follows, that the starburst 
luminosity would have to be ~ 100 times larger than the 
maximum starburst luminosity Lm given by MQT if the 
onset of a photon-driven wind is to limit the luminosity of 
a galaxy. If this were the case, the bolometric luminosi- 
ties of galaxies would be ~ 10 — 100 times larger than 
their respective quasars, which clearly violates current 
observations. 

More recently, Thompson et al. (2005) postulate that, 
on the basis of observations of bright local starbursters 
such as Arp 220, prompt starburst photons self-regulate 
the star formation rate via the infrared dust opacity. If 
the dust is optically thick to far IR emission, this model 
has many of the benefits of an optically thick CR feed- 
back model. Specifically, the net force imparted by each 
photon is amplified due to multiple scatterings and the 
momentum imparted is insensitive to the source distri- 
bution. In order to be optically thick to the far IR, large 
galactic gas surface densities (£ s > 1) and large dust 
temperatures (T^ ~ 100 K) over the bulk of the star- 
burst region, as indicated by some observations of Arp 
220, arc required (Downcs & Solomon 1998; Soifer et al. 
1998). In this case, it is possible that the infrared pho- 
ton pressure may exceed the CR proton pressure near 
the midplane, especially if CR destruction becomes im- 
portant (see §3.4p . However, this interpretation that the 
nuclear activity in Arp 220 is powered by star formation 
has recently been called into question (e.g. Downes & 
Eckart 2007). 

2.4. Aspects of Galactic Scale CR-Driven Winds 

As stated in t )2.1l starbursters display outflows with 
velocities significantly in excess of the escape speed v esc 
or the velocity dispersion a. Interestingly, line-driven 
winds of massive stars and dusty radiation-driven winds 
of asymptotic giant branch stars also display fast winds 
with velocities in excess of v csc due to multiple scatter- 
ings, outside the main body of the respective object. 

Radiation (or cosmic radiation) pressure driven winds 
occur when the luminosity of the source approaches and 
surpasses the frequency-averaged Eddington limit i.e., 
when the radiation force exceeds gravity. If the radia- 
tively supported source is in quasi-hydrostatic balance 



beneath its photosphere, then a wind can only develop if 
the opacity increases above the main body's photosphere. 
For example, in the winds of AGB stars, the fact that 
dust sublimates at temperatures below the photospheric 
temperature leads to an opacity increase with decreasing 
density (Salpeter 1974; Goldrcich & Scovillc 1976; Ivezic 
& Elitzur 1995). 
To see how this works, we write 

kF F dv 

P = 7-=P.9 + P«-r ( 23 ) 

c \c dr 

where k, A, F, v, and g is the radiation opacity, mean free 
path, flux, outflow velocity, and gravitational accelera- 
tion, respectively. If k increases with decreasing density, 
then the ratio of radiative to gravitational acceleration 
increases as a parcel of gas moves outwards. Integrating 
over volume allows us to write 

/oo / p \ L 

r 2 dr [j- -pgj ~ -t w = M Voo (24) 

i.e., infinite acceleration may occur in the limit of infi- 
nite optical depth. Here, t w is the optical depth of the 
wind, M = 4:TrR 2 pv is mass loss rate, Voo is the terminal 
velocity, and Rq is the launching radius. In the above 
expression, the gravitational acceleration is neglected in 
the second term since it is assumed to be small in com- 
parison to the radiative acceleration. Of course, t w is 
limited by energy conservation, which we discuss in Ap- 
pendix ini 

The scenario outlined above for radiatively-driven 
winds of stars may analogously occur for CR-drivcn 
galactic winds. Cosmic ray transport in the Milky Way 
is often modeled with a constant diffusivity over a large 
rarefied kpc halo; a picture that is consistent with obser- 
vations of the galactic CR distribution . Since the diffu- 
sivity varies little - or not at all - with density, the CR 
cross section and thus opacity, increases with decreasing 
density. 

In §3] we discuss how the CR mean free path and cross 
section directly depend only upon the magnetic field fluc- 
tuations at the Larmor scale. That is, the level of cou- 
pling between CRs and matter is, in principle, indepen- 
dent of column density - an ideal situation for winds 
with multiple scatterings. For example, if the streaming 
instability is the source of resonant magnetic fluctuations 
(see ^3.3.ip . then kcr oc p~ 1//2 /B, where kcr is the CR 
"opacity." 

As stated in §2.3.21 UV photons released from the 
surface of massive stars are absorbed and then quickly 
downgraded into effectively non-interacting IR photons. 
Therefore, t 7 ~ 1 for photon feedback - if one accepts an 
idealized one-zone approximation (see §2.3.2| . When the 
momentum of the radiation field is so weakly coupled to 
the matter, it is difficult to imagine how photon-driven 
feedback leads to wind velocities in excess of v csc by a 
factor of ~ a few, unless the light to mass ratio exceeds 
that of eq. ([18]) by a factor of ~ 10 or so, since the max- 
imum luminosity Lm on v 2 ^ for a fixed enclosed mass M 
according to formalism of MQT (see e.g., their eq. (26)). 

Independent of galactic mass, SNe shocks heat inter- 
stellar gas to temperatures of about ~ IkeV, which cor- 
responds to a sound speed of ~ 300 km/s (Martin 1999, 
2004). If SNe lead to a thermally-driven wind that is re- 
sponsible for ejecting the interstellar medium, then the 
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outflow velocity is comparable the sound speed. It fol- 
lows that thermally-driven winds from SNe are unlikely 
candidate for explaining fast > 300 km/s galactic winds. 

2.4.1. Upper Limit on r CR from Momentum and Energy 
Conservation 

Even in our Galaxy, the interpretation of direct and 
indirect measurements of CR transport may lead to qual- 
itatively different pictures as to how CR diffuse through 
the Milky Way and its corona-like halo. Nevertheless, we 
may set some rough limits on the CR optical depth t cr 
without knowing the details of the particular mechanism 
responsible for the resonant scatterings. 

In the Appendix |Dj we briefly estimate an upper limit 
on t cr , following a well established result from the stel- 
lar winds literature (Lamcrs & Cassinelli 1999; Ivezic 
& Elitzur 1995; Salpeter 1974). For optically thick 
radiation-driven winds, momentum and energy conserva- 
tion set an upper limit for the wind optical depth roughly 
as r CR ~ c/v w , where v w is the wind terminal velocity. 

The maximum value of t cr takes on interesting val- 
ues when v w is large enough as to merit escape from 
a given galaxy's gravitational potential. As previously 
mentioned, since CRs inject roughly ~ IkeV/baryon into 
the interstellar medium, the deepest gravitational poten- 
tial from which they can unbind gas corresponds to a 
velocity dispersion a ~ 300 km s -1 . If v w ~ a, then the 
corresponding optical depth is t cr ~ 10 3 , which is close 
to the Milky value. 

3. PHYSICAL PROPERTIES OF INTERSTELLAR COSMIC 

RAYS 

For lack of a better choice, we use the Galaxy in or- 
der to calibrate the properties of cosmic ray injection, 
transport, and star formation in general. 

Cosmic ray protons within the relatively modest en- 
ergy scale of ~ 1 — 10 GeV arc responsible for ~ 90% 
of the Galactic CR pressure near the solar neighbor- 
hood. Above 10 GeV, The CR energy energy spectrum 
n(E) is well described by a single power law such that 
n(E) cx E~ n from - 10 - 10 6 GeV. As already implied, 
the CR spectrum is steep such that n > 2, with a com- 
monly quoted value of n ~ 2.7. 

The spatial distribution of Galactic CRs is highly 
isotropic such that the energy-weighted CR mean free 
path A CR ~ 0.1 — 1 pc. This remarkably low value for 
A CR is most likely the result of CRs scattering off of inter- 
stellar magnetic irregularities. Then, for typical Galac- 
tic scale height H ~ 1 kpc, the interstellar matter of the 
Galaxy may be considered as possessing an optical depth 
to CRs r CR ~ 10 3 — 10 4 . Therefore, the CR pressure P CR 
is related to the CR flux F CR by the simple relation 

P CR ^F CR ^f (25) 

over a CR pressure scale height. 

In the Galaxy, a core-collapse supernova occurs once 
every century, leading to a core-collapse supernova lumi- 
nosity L S N 

L SN ~ 10 51 erg/10 2 yre - 3 x 10 41 m SF ergs -1 (26) 

where m SF is the normalized star formation rate in units 
of IMq/yt. The injection of CRs results primarily from 
the interaction of supernova shocks with the interstellar 



medium (Ginzburg & Syrovatskii 1964). However, CR 
deposition from other types of starburst outflows, such 
as stellar winds from massive stars, may also significantly 
contribute (Schlickeiscr 2002). In terms of the star for- 
mation rate m SF and an CR efficiency e CR , we may write 
the CR luminosity as 

^cr = £ cr™sf c2 = 6 x !0 4 ° £ 6 rh SF ergs -1 . (27) 

where e$ is equal to e CR in units of 10 -6 or in other words, 
L CR is normalized to ~ 20% of L$n- 

3.1. Observational Constraints on A CR from the Leaky 
Box and Galactic Halo Models 

Upper limits on the anisotropy of the CR distribu- 
tion function, the relative abundances of primary to 
secondary CRs, and the abundance ratios between the 
radioactive isotopes of a given CR species sum up 
the directly available observations of CR proton trans- 
port (Schlickeiscr 2002). Indirect tracers of CR be- 
havior include radio emission, presumably powered by 
brcmsstrahlung and synchrotron losses of CR electrons, 
as well as gamma-ray emission resulting from the decay 
of mesons produced in proton-proton collisions. Alto- 
gether, these observational constraints indicate that CRs 
stay in contact with the Milky way for time scales in the 
neighborhood of 10 7 yrs. 

Though the direct and indirect observations of CR 
transport dictate that CRs scatter t? r ~ 10 6 — 10 8 times 
before leaving the Milky Way, they do not uniquely de- 
termine the volume of the region in which they are con- 
fined (see e.g., Ginzburg et al. 1980 for a discussion). 
In the one-zone model, or "grammagc" formulation, two 
important pieces of information can be extracted; the 
galactic residence time for a CR proton of a given energy 
and the average hydrogen number density (hh) during its 
encounters with the interstellar medium. For ~ 1 — 10 
GeV CR protons (tih) ~ 0.2 cm -3 - significantly smaller 
than average value for the ~ 100 pc scale height gaseous 
disk, where is tih ~ lem . In the grammage formula- 
tion, typical residence timcscalcs in the disk are of order 
~ 3 x 10 6 yrs. Together, this implies an optical depth 
in the disk of r CR ~ 10 4 and a mean free path of order 
A CR ~ O.lpc. 

The average hydrogen number density hh ~ 1 cm for 
the 100 pc disk, which implies CR protons must somehow 
avoid regions of relatively large density while maintain- 
ing a small level of spatial anisotropy, of order one part 
in 10 3 - 10 4 . The mass distribution of the Milky Way 
is partitioned into several phases of various densities and 
temperatures. Even if CR protons mainly restrict them- 
selves to a phase with characteristic hydrogen number 
density njj ~ 0.2 cm -3 , it seems unlikely that such a 
phase of the interstellar medium is homogeneously dis- 
tributed in space at the one part in 10 3 level, in the solar 
neighborhood. 

In light of the difficulties with the one-zone interpreta- 
tion of grammagc measurements, another scenario that 
is often used to interpret grammage measurements is 
the disk + halo model (Ginzburg & Syrovatskii 1964; 
Ginzburg et al. 1980). By modeling the vertical pro- 
file of the Galaxy, the quantitative interpretation of a 
CR proton's trajectory is much different in comparison 
to one-zone disk models. Here, 1-10 GeV CR protons 
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are thought to be rather uniformly distributed in a dif- 
fuse ~ 1 kpc disk with typical hydrogen number densities 
njj ~ 0.1 cm~ 3 , with a characteristic escape time of or- 
der ~ 10 7 yrs (Strong & Moskalenko 1998; Webber & 
Soutoul 1998) . In this picture the CR mean free path 
^CR ~ 1 P c - 

3.2. Constraints on ecR and A CR from other Galaxies 

Although the mechanism responsible for the 1-10 GeV 
cosmic ray mean free path in the Milky Way is not com- 
pletely understood, it seems that the value of the cosmic 
ray mean free path Acr may be roughly the same from 
galaxy to galaxy, over a large range of star formation 
rate. Important constraints on cosmic ray propagation 
come from studies of the tight far-infrared/radio corre- 
lation, which spans over many decades in star formation 
rate (see e.g. Condon 1992). The far- infrared (FIR) 
emission is due to optical/UV star burst photons being 
reprocessed and down-graded into the FIR by dust, while 
the GHz synchrotron emission is powered by 1-10 GeV 
cosmic ray electrons. Like cosmic ray protons, cosmic ray 
electrons are supposedly accelerated at the shock fronts 
of supernovae and stellar winds, by a common Fermi ac- 
celeration mechanism. 

Above an energy of 1 GeV, where both electrons and 
protons are relativistic, the CR mean free path of both 
species are the same as long as the scattering opacity 
results from small scale resonant magnetic irregularities. 
Therefore, a measurement of the 1-10 GeV CR electron 
mean free path serves as a good proxy for the 1-10 GeV 
CR proton mean free path. By setting the value of the 
magnetic energy density approximately equal to the to- 
tal cosmic ray energy density - the so-called minimum 
energy argument - the total CR electron energy content 
within the synchrotron cooling radius is extracted. Along 
with the total synchrotron power within the cooling ra- 
dius, this yields the synchrotron cooling time for the 1-10 
GeV CR electrons. 5 Together with the cooling time, the 
characteristic emission scale allows for a determination 
of the the 1-10 GeV CR mean free path, which is typi- 
cally close to 1 pc or in other words, close to Milky Way 
value (Condon 1992). 

As previously stated, the interaction of energetic out- 
flows from massive stars (supernova and stellar winds) 
with their ambient gaseous environment is commonly 
thought to be the site of interstellar CR production. In 
our galaxy, roughly 10% of the core-collapse supernova 
power is converted into 1-10 GeV CR proton power, £ CH , 
and of that only ~ 1% of L CR manifests itself in 1-10 
GeV CR electron power. Under the assumption that all 
of the GHz radio emission from starforming galaxies re- 
sults from synchrotron cooling of CR electrons and that 
all of the FIR emission is due to reprocessed starlight 
onto dust, the tight FIR-radio correlation indicates that 
the efficiency of producing CR electrons, per unit of star- 
formation, is remarkably constant from galaxy to galaxy. 

5 We do not understand why the total (electron+proton) cosmic 
ray energy density is utilized in the traditional minimum energy 
argument, since synchrotron emission is produced only by cosmic 
ray electrons. Nevertheless, the assumption that the magnetic en- 
ergy density scales with the cosmic ray proton energy density does 
not appreciably alter the value of the synchrotron cooling time as 
long as the the ratio of cosmic ray proton to electrons does not 
significantly exceed the Galactic value. 



Finally, if a magnetized Fermi shock acceleration process 
is responsible for producing interstellar CRs from a bath 
of thermal particles at a common sub-relativistic energy, 
then one expects on theoretical grounds that the ratio 
of CR protons to electrons to be ~ 100 or so, close to 
the observed value in the interstellar medium (Schlick- 
heiser 2001). Altogether, this implies that the efficiency 
of producing 1-10 GeV CR protons e CR varies little from 
galaxy to galaxy. 

3.3. The Source of Small Scale Resonant Magnetic 
Irregularities or: A Great Uncertainly in Cosmic 
Ray Physics 

3.3.1. Streaming Instability 

CRs flow freely along large scale laminar magnetic 
fields as they tightly gyrate about their guiding centers. 
In the presence of small scale magnetic fluctuations, the 
CRs themselves can amplify waves in which they are res- 
onant, a phenomena known as the streaming instabil- 
ity (Kulsrud & Pearce 1969; Wentzel 1974). The phase 
velocity of the waves in question (e.g. fast or Alfvenic 
waves) is much smaller than the speed of light, the ve- 
locity at which the CRs propagate. Thus, a CR perceives 
a wave packet of magnetic disturbances as stationary in 
time, but not space (Kulsrud 2005). Therefore the ap- 
propriate resonance condition for the interaction between 
a CR beam centered in energy at some Lorentz factor 
7 and a wave packet of magnetic disturbances of some 
width Ak is that the Larmor radius of the CRs equal 
to the component of the wavelength projected along the 
laminar field. 

The growth time of the streaming instability is ex- 
tremely short in comparison to local dynamical times 
of interstellar matter for 1-10 GeV CRs. In the absence 
of damping CRs excite waves which propagate along the 
background field with a growth rate proportional to the 
ratio of the streaming or diffusion velocity to the Alfvcn 
speed vd/va- As the amplitude of the waves grow, for- 
ward momentum is extracted from the streaming CRs 
as they are isotropized, while simultaneously reducing 
their streaming velocity. Wave excitation continues un- 
til vd ~ va, where the instability is quenched and the 
instability's source of free energy has been removed. 

Now, we may write the CR optical depth as 



From this, the CR optical depth of the Milky Way halo 
with va ~ 150 km s^ 1 is r CR ~2x 10 3 , which is roughly 
the observed value. 

In order for the streaming instability to operate at 
given resonant wavelength, the growth rate of the in- 
stability must surpass the individual damping rates re- 
sulting from all sources of dissipation, such as ion-neutral 
drag. The growth rate due to the streaming instability 
for a given resonant Alfven wave is proportional to the 
number density of CRs n[E) resonant with the magnetic 
disturbance in question. Due to the steepness of the CR 
energy spectrum, the CR number density at 100 GeV 
n(E = 100 GeV) is sufficiently low that driving resulting 
from the streaming of CRs cannot overcome ion-neutral 
drag (Cesarky 1980). The argument presented above rep- 
resents a major shortcoming of the "self-confinement" 
picture of CR transport. 
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3.3.2. Kolmogorov Cascade 

On its own, the interstellar gas may be a rich source of 
small scale resonant Alfvenic disturbances. In this case, 
the CR mean free path A CR is limited by the amplitude 
of magnetic fluctuations on the Larmor scale. 

In the frame of the CR beam, magnetic irregulari- 
ties of the interstellar medium are effectively stationary 
and thus, diffusion occurs only in space since energy is 
roughly conserved per scattering event. Spatial diffusion 
proceeds in pitch-angle with disturbances that are reso- 
nant with the individual CRs in momentum-space. The 
appropriate diffusion coefficient T>qr is roughly given by 



V C r ~cr L 



D- 



SB 2 (r L )' 



(29) 



Here, SB (r^) is the amplitude of magnetic fluctuations 
on the Larmor scale. If the magnetic fluctuations follow 

1/3 

a Kolmogorov scaling, then SB = SB (A/A ) , where 
A, Ao, and SBq is the eddy scale, the stirring scale, and 
the magnetic fluctuation amplitude at the stirring scale. 

The above expression reflects the notion that a given 
CR can only wander in pitch-angle of its guiding center 
with the aid of resonant magnetic irregularities. Since 
the phase of the irregularities are taken as random, a 
CR must execute B 2 jSB 2 {ri,) gyrations before its net 
trajectory, and thus momentum, is altered at the or- 
der unity level. It follows that the CR mean free path 
A CR ~ 2?cr/c, when the source of magnetic deforma- 
tions is due to turbulence. With this, the CR optical 
depth in the presence of resonant Kolmogorov magnetic 
turbulence reads 



#Mcr ' 



cH/V CR 

1/3 / jj\ 2/3 
V 



(30) 



where £b = SB /B. 

As previously mentioned, the Alfvenic cascade must 
possess, with sufficient amplitude, eddies that are reso- 
nant with the CRs at the Larmor scale. In order to be 
resonant, the Alfven waves should propagate primarily 
along the field. However, it is now generally accepted 
that, in accordance with Goldrcich, Sridhar and Lith- 
wick's (Goldreich & Sridhar 1995; Lithwick & Goldreich 
2001) theory of interstellar MHD turbulence, Alfvenic 
power flows only into eddies that propagate in the direc- 
tion perpendicular to the local large scale field. Thus, for 
Ao/^l 3> f , there may be virtually no resonant turbulent 
power left at the Larmor scale. 

3.4. Cosmic Ray Destruction 

Cosmic rays can be destroyed once they interact with 
an ambient proton, subsequently producing pions. At a 
few GeV, the cross section for this inelastic process is 
~ 30 mb. Therefore, the characteristic destruction time 
tn for GeV CR protons is (Schlickeiser 2002) 



108 

n H 



-yrs, 



(31) 



which should be compared to the diffusion (or escape) 
time over a galactic scale height H that is occupied with 
gas of average density riu- 



In the Milky Way, the destruction timescale t n does not 
play an appreciable role in the transport of CR protons. 
Note that this statement is independent on whether or 
not the one-zone or disk + halo model is utilized to inter- 
pret grammage measurements since the hydrogen num- 
ber density throughout the Milky Way is relatively low. 

In young bright starbursting galaxies may reach 
values of 10 3 — 10 4 cm~ 3 , corresponding to a destruc- 
tion time of tir ~ 10 5 — 10 4 yrs. It follows that in some 
of the interesting cases in which CR feedback might be 
important, CRs may in fact be destroyed before they tra- 
verse a galactic scale height. In this case, it is useful to 
employ an effective optical depth when diagnosing the 
level of feedback that CRs exert back upon the inter- 
stellar medium. If a CR is destroyed before it escapes 
the galaxy, it still has the the opportunity to exchange 
momentum with a given galaxy's gaseous component. 

In order to understand the effects of CR destruction, 
we have calculated the vertical variation of CR pressure 
for a homogeneous disk geometry in Appendix [Cj As in 
the case of scattering dominated, radiative atmospheres 
(see e.g. Felten & Rees 1972), the physics is easily un- 
derstood in terms of random walk arguments. We define 
the effective CR optical depth, r CR , as 



tdiff 



(32) 



where r CR ~ H/X CR , t^ig is the diffusion time for the CR 
protons over the galactic scale height H, and is CR 
the "absorption" optical depth due to pion production. 

The above expression acknowledges the fact that CRs 
must scatter t 2 r times in order to traverse a distance 
H in a diffusion time idiff- In the limit that t„ < t^g, 
the square of the effective optical depth T 2 is equal to 
the number of times a CR proton scatters and exchanges 
momentum with the galaxy before its destruction. 

For sufficiently dense and extended systems, destruc- 
tion of CR protons due to pion production will inhibit 
the efficiency of CR feedback. Consider the example 
of a gas-rich disk with n# ~ 10 4 cm" 3 and a disk 
scale height H ~ 100 pc. The destruction timescale 
t v ~ 10 4 yrs is quite short. If r CR ~ 10 3 for this 100 
pc disk then A CR ~ 0.1 pc and idiff ~3x 10 5 yrs. There- 
fore, t cr - 0.2 r CR - 200. 

For galactic environments in which 1-10 GeV CRs are 
destroyed by hadronic processes before their escape eq. 
(|3"2"|) provides a prescription for estimating the CR con- 
tribution to hydrostatic balance. If we assume the CR 
source is independent of height in the disk, then pressure 
gradient (force) at the surface is given by 

dP Ft F n 

u± CR ± CR ' CR ± CR / (33) 

dz He A CR c V tdis 

Deeper within the atmosphere the force is smaller due to 
cancellation of momentum between upward and down- 
ward moving CRs, going to zero (by symmetry) at the 
midplane. 

In situations were the CR mean free path A CR is small 
and the hydrogen number density nn is high, some 1- 
10 GeV CRs may never escape the galaxy since meson 
production will destroy them. The typical vertical dis- 
placement a CR drifts before its demise ^ CR is given by 
l li 

l CR ~ A CR iV s , where N s = t n /t s is the number of scat- 
tering before destruction and t s = A CR jc is the time in 



Cosmic Ray Feedback 



9 



between scattering events. Rewritten in terms of Z CR , eq 
((33|) becomes 



dz 



A, 



(34) 



where we have implicitly assumed l CR < H (i.e. CRs 
are typically destroyed before leaving the galaxy) . In the 
absence of CR destruction processes, the term in paren- 
theses would be unity. The effect of CR destruction has a 
simple interpretation. Only those CRs produced within 
l CR reach the surface before being destroyed. Therefore, 
only the fraction l CR /H of the CRs contribute to the 
force at the surface. 

A useful estimate of l CR /H can be obtained from the 
relation tih ~ Y> g /(rn p H) where T, g is the gas surface 
density. We then find that 



CR 

H 



60 g cm 2 A c 



H 



1/2 



(35) 



It is clear from this expression that if either S g is suf- 
ficiently large and/or \cn/H is sufficiently small, most 
CRs will be destroyed before escaping the star burst. Of 
course, reducing Acr at fixed H produces a net increase 
in the force due to the dependence in eqs. ([33|) and 
(|34|) . However, this force is larger only in an increasingly 
smaller fraction (l CR /H) of the volume near the surface. 
Deeper in the atmosphere, a CR feedback scenario will 
run into some of the same problems photon-driven feed- 
back suffers from as mentioned in Appendix [Bj Near the 
surface there will be a force comparable to our estimate 
in eq. (|33|) . but throughout most of the starburst vol- 
ume the force due to upward moving CRs from any one 
source is mostly canceled by downward moving CRs from 
neighboring sources above. 

Finally, we note that in their eq. (44), Jubelgas et al. 
(2006) include the effects of CR feedback by assuming a 
pressure pressure P CR cx t„ . This gives a sharp reduction 
in pressure with increasing density since t v cx . As 
discussed in Appendix [Cl this assumption is appropriate 
deep within galaxy if i CR is much less than the charac- 
teristic galactic radius. However, in the outer regions 
(within ~ Z CR ) the CR forces will be significantly larger 
due to the dP OH j dz cx ^/t^ dependence in eq. (|33|) . 

3.5. Overall Picture and Summary 

If in fact the expulsion of CRs is the mechanism respon- 
sible for limiting the luminosity of star bursting galaxies, 
then the CRs must be well-coupled to the interstellar 
gas such that the CR optical depth r CR is large. Of 
course, there are no direct measurements of CR proper- 
ties in other galaxies, which is why we have primarily 
utilized observations and theories of the CR transport in 
the Milky Way as our benchmark. Even for the Milky 
Way, the value of the CR mean free path A CR extracted 
from the data may vary by an order of magnitude, where 
A CR ~ Ipc if the halo model is utilized and A CR ~ O.lpc 
for the leaky box interpretation. 

From the theoretical side, the two most commonly ac- 
cepted explanations for the source of small scale resonant 
magnetic perturbations possess severe weaknesses. Since 
the growth rate of the streaming instability at a given CR 
energy is proportional to n{E) cx E~ 2J ', the streaming 
instability is not strong enough to overcome ion-neutral 



damping at energies above 100 GeV (Cesarsky 1980). 
Perhaps one might account for the lack of an-isotropy - 
and thus small mean free path - above 100 GeV asserting 
that the source or magnetic irregularities on those scales 
comes about from a turbulent cascade of resonant mag- 
netic power. However, it seems unlikely that the slope of 
the CR distribution function with energy remains con- 
stant over ~ 6 decades in energy if the streaming insta- 
bility were to inject power at energies below ~ 100 GeV. 

Despite these disturbing inconsistencies in the overall 
picture of CR transport in the Milky Way, what we know 
for certain is that the scattering optical depth r CR is large 
such that t cr ~ 10 3 — 10 4 . Furthermore, we may con- 
clude, with an acceptable degree of confidence, that the 
origin of the small CR mean free path A CR ~ 0.1— lpc in 
the Milky Way is magnetic in nature. Within the context 
of galaxy formation, the fundamental issue is constrain- 
ing how r CR and A CR evolves with environment (redshift 
and halo mass). In what follows, we briefly discuss the 
manner in which r CR might evolve from galaxy to galaxy. 

3.5.1. Does Cosmic Ray Feedback Result in the 
Faber- Jackson Relation? Maybe. 

Eq. (|13|) appears strikingly similar to Faber- Jackson 
(1976) relation. Under the plausible assumption that 
galaxies residing on the Faber-Jackson relation are faded 
and elderly versions of their bright former starbursting 
past, then a momentum-driven feedback mechanism may 
indeed lead to the Faber-Jackson relation. If CR feed- 
back is responsible for the relation, then the CR optical 
depth t cr can only evolve weakly with velocity dispersion 
a. 

The CR optical depth r CR depends on galactic mag- 
netic field strength if either the streaming instability or 
a Kolmogorov cascade is responsible for the required res- 
onant magnetic fluctuations. The situation is more com- 
plicated in the case where a magnetic cascade is respon- 
sible for the CR optical depth since knowledge of £b , Ao , 
H and the power law index of the cascade (= 1/3 for 
Kolmogorov) is required as well. At the moment, a first 
principles calculation of r CR in the Milky Way would be 
a difficult task for reasons already mentioned. Certainly, 
an attempt to determine t cr from basic physical argu- 
ments in other galaxies is out of reach. 

As stated in §2.4.11 t cr cannot assume and indefinitely 
high value. In other words, only a finite number num- 
ber of scatterings are required before the CR beam van- 
ishes due to adiabatic losses and thus, energy conserva- 
tion enforces an upper limit on r CR . Furthermore, if the 
CR residence time, and consequently r CR , becomes too 
large then meson production will limit the effective opti- 
cal depth t cr to a moderate value. Note that the r CR is 
bounded from above as a result of constraints that have 
nothing to do with the ultimate - and highly uncertain - 
mechanism that generates the required resonant Alfvcnic 
disturbances that scatter the CRs. 

On less solid grounds, we can conclude that r CR cannot 
take on values that are too low for a given galaxy. If 
the galactic magnetic field is too weak, then the growth 
rate of the streaming instability is enormous since the 
ratio c/va becomes very large. On the other hand, if the 
galactic field is quite strong on the stirring scale, then 
one would expect that the amount of resonant power at 
ri, is quite large, leading to a relatively small mean free 
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path. 



4. CONCLUSION 



From the process of star formation, cosmic rays rep- 
resent only a minor form of energy release. Yet, their 
momentum coupling with the interstellar medium is ab- 
solute as result of their small scattering length. As long 
as the CR mean free path A CR does not greatly differ 
from the inferred value in the Milky Way, the collec- 
tive force that CRs exert upon the interstellar medium 
as they leak out is great enough to support the gaseous 
component against gravity and may even unbind it. 

If CR feedback is an important mechanism in deter- 
mining the maximum luminosity of a starbursting galaxy, 
then the CR optical r CR must be large in order to com- 
pensate for the relatively small CR luminosity. Con- 
sequently, galactic CRs necessarily impart momentum 
upon the interstellar medium in the direction opposite 
to the galaxy's gravitational center. Thus, CR feedback 
is highly independent of the spatial distribution of CR 
injection sites such as SNe and stellar winds. 

In analogy with the photon-driven winds of massive 
and asymptotic giant branch stars, CR-driven galactic 
winds may reach terminal velocities in excess of the es- 



cape velocity by a factor of ~ a few. Due to multiple 
scatterings in the galactic halo, above the main body of 
the galaxy, the CR beam may indefinitely accelerate a 
parcel of fluid up to the point where all of the CRs are 
rcdshiftcd away. 

The main uncertainty in the CR feedback theory pre- 
sented here is in determining the manner in which the CR 
mean free path A CR (or optical depth t cr ) and radiative 
efficiency e CR varies from galaxy to galaxy. The optical 
depth r CR cannot approach arbitrarily high values since 
adiabatic losses and meson production will eventually re- 
move energy out of the CR beam. 
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APPENDIX 

APPENDIX A: A GALAXY AS AN ISOTHERMAL SPHERE 

As a baseline model, we assume the mass profile of a given galaxy is that of a singular isothermal sphere. This 
assumption is more accurate for hot stellar systems such as elliptical galaxies and bulges and less accurate for disk 
galaxies. The density profile of the gaseous component is given by 



which corresponds to an enclosed gas mass of 



p(r) 



M g (r) 



i> 2 

2nGr 2 



G ' 



(Al) 



(A2) 



In the above expressions, f g , a, and r are the gas fraction, velocity dispersion, and galactic radius, respectively. 

APPENDIX B: BREAKDOWN OF THE ONE-ZONE MODEL FOR RADIATIVE FEEDBACK 

We assume that every massive star is embedded such that every single one of its emitted UV photons are absorbed 
by dust grains and reprocessed into infrared radiation close to the star. This is not the case in the Milky Way where 
massive stars are embedded for only ~ 10% of their lifetimes (Wood & Churchwell 1989). However, for starbursting 
sources, the ratio of infrared to UV radiation is large, implying that our assertion of every massive star being embedded 
is correct. 

UV photon power is converted into the IR at a characteristic radius of ~ 0.1 — lpc from the center of each individual 
star. It is in this region, the so-called "cocoon," in which momentum is directly exchanged between light and matter. 
Note that the characteristic hydrogen number densities in this region range from n ~ 10 5 — 10 7 cm -3 (Koo & Kim 
2001; Churchwell 2002). 

There are ~ 10 4 O stars in the Milky Way. For a starbust galaxy with m SF ~ 100 or so times the Milky way value, 
one expects there to be No* ~ 10 6 or so massive stars responsible for both the starburst luminosity as well as the 
injection of radiation momentum. If the volume of this galaxy is ~ 10 9 pc 3 , then the typical separation between each 

1/3 

O star is Ax ~ (V/No*) ~ 10 pc for a homogeneous and spherically symmetric distribution of O stars. 

Furthermore, the photons from each O star may be thought of as being responsible for injecting momentum into 
l/A"o* of the given galaxy's gaseous component. It follows, that momentum balance for a parcel of dusty gas placed 
somewhere between two typical massive stars may be expressed as 



o* 



47rAa 



Ax/2' 



(Bl) 



where kuv and Lo* is UV opacity on dust and the characteristic luminosity of an O star. In the above expression, 
the radiation force imparted at the surface of the dusty "cocoon" is matched by a non-linear momentum sink, which 



Cosmic Ray Feedback 



11 



can be thought of as either an outward moving shell or a parsec-scale turbulent eddy. Integrating over the volume in 
between the two stars gives 

Ln* v 2 

f T ~ M ^ (B2) 

where Ma ~ f g M(R)/No* is the fraction of gaseous galactic mass in which a given massive star is responsible for 
accelerating. Noting that Lo* = ^sf/^Vo* together with the "single scattering approximation" (r = 1) of MQT, we 
may write 

2 L SF G (Ax/2 




*~iiu*\rty (B3) 

For Lsf = Lm, where Lm = 4/ 9 ca 4 /G is maximum stellar luminosity derived in MQT, we find that 

' (B4) 

For a bright starbust galaxy with No* ~ 10 6 , the radiation force may lead to random fluid motions in the gas with 
characteristic velocities of only ~ 1/10 that of the isothermal velocity dispersion. Therefore, it seems highly unlikely 
that radiation-driven feedback can lead to the unbinding of a galaxy's gaseous component. The clustering of massive 
stars won't alleviate this major concern. Even if O stars are embedded together 100 at a time behind a common dust 
sublimation layer, the random velocity v increases by only a factor of ~ 2. 

COSMIC RAY TRANSFER WITH DESTRUCTION: A RANDOM WALK MODEL FOR THE COSMIC 

RAY PRESSURE GRADIENT FORCE 

We begin with the equations relating CR pressure and flux for a plane-parallel disk-like geometry. First we assume 
that CR flux Fen is determined by 

dF 

^L = -i-i 3 P CR + Q CR . (CI) 

Here we have assumed 7cr = 4/3, t„ is the CR lifetime to pion production (proton-proton scattering), P CR is CR 
pressure, and Q CR is the CR source term. Other processes may destroy CRs, but we expect pion production to be 
most important for starforming galaxies. In the diffusion dominated limit F CR has a simple relation to P CR 

dP F 

— ^ = — -2*^. (C2) 
dz c A CR 

Combining cqs. IjCip and (|C2[) . and using the definition X v = t v c we find 

d 2 P 3 P (1 

CR OJ CR ^CR 



dz 2 A cr Att A 



CH ' 



This equation now resembles the angle averaged radiative transfer equation assuming the Eddington relation P CR = 
E CR /3 applies. Following the normal prescription for a scattering dominated radiative atmosphere (see e.g. Rybicki 
& Lightman 1979), we now define an effective mean free path l CR = (A CR A^/3) 1 / 2 and <fr* = — l~ R dz. This yields the 
familiar form 

For simplicity we will assume X^Q^ is constant. With this assumption Q CR = F max /H where H is the disk scale 
height and F max — L OK /(4nR 2 ). The solution to eq. (|C3|) is then given by 

In order to specify P+ and P_, we need two boundary conditions. First, we will assume that disk has a midplanc 
defined by r m = H/X CR where P CR vanishes due to symmetry. At the surface, we assume there is no incoming CRs. 
In the two-stream approximation (see e.g. Rybicki & Lightman 1979) this gives 

1 dPc* _ p 



\/3 dr CR 

at r CR = 0. With these assumptions, it is easy to show that 

exp(—T*) + exp(r* — 2H/l OK ) 



p _ AttQqr 

CR ~ 3c 



1 - 



l + ^+exp(-2H/l CR )(l-^) 



(C4) 
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Fig. CI. — Cosmic ray pressure as a function of atmospheric depth for varying levels of cosmic ray destruction. The midplane corresponds 
to z = 0. The vertical axis is the cosmic ray pressure in units of H F msx /c\ CR . The ratio of the cosmic X CR /H = 0.1. The solid line 
corresponds to pure scattering and no destruction i.e., H/l CR = 0, while the H/l CR = 0.5 (dotted), 1.0 (dashed), 3.0 (dot-dashed) and 5.0 
(triple-dot dashed). For relatively large levels of CR destruction, the midplane pressure is reduced by a factor ~ l^ R /H 2 , as in the case of 
Jubelgas et al. (2006). 



where we have defined e = A CR /A 7r (not to be confused with supernovae efficiency) for convenience. Note that, by 
assumption, f«l. 

A number of results are immediately apparent from eq. (|C4[) . First, we confirm that P CR is a maximum at the 
midplane. Furthermore, in the limit H ^> l CR where P CR = t 7r Q cn /3. This simply shows that if all CR arc destroyed 
before they escape the total energy density of CRs is equal to the rate at which they are produced multiplied by the 
time it takes to destroy them. Note that this is the assumption of Jubelgas et al. (2006). 

Ultimately, we are interested in the force (per unit volume) the CRs provide, which is given by 

dP C n = ^tQcr exp(-T*)-exp(T*-2H/l CR ) 
dz 3cl CR 1 + y/e + exp(-2H/l CR )(l- y/e)' 1 ' 

It is clear from the form of the numerator that this force vanishes at the midplane. This must be the case, since the 
CR flux vanishes at the midplane due to the assumed symmetry. 

Now consider the force near the surface. We'll first examine cq. (|C5[) in the limit that H l CR . In this case the 
numerator and denominator simplify considerably and we find 

dz ~ 3d L z ' 

For this case we see that the force increases linearly with height z above the midplane. The maximum force is reached 
at the surface z ~ H where 



dz 3ll R c A CR c A CR ( 




Fig. C2. — Cosmic ray pressure gradient as a function of atmo sphe ric depth for varying levels of cosmic ray destruction. The values for 
all of the input parameters are identical to those chosen in FieurclCTI The pressure gradient and therefore, the vertical cosmic ray pressure 
gradient force vanishes at the midplane in all cases, due to symmetry. Furthermore, the cosmic ray force and pressure gradient attain 
their maximum values at the top of the atmosphere, whether or not destruction is important. For the cases in which CR destruction is 
important, the maximum value for the CR force is reduced by l CR /H. It follows that the Eddington limit in cosmic rays L Bdd CR increases 
by a factor of H/l CR relative to the pure scattering case, for situation where CR destruction is important. 



which is the expected result if CR destruction is negligible. 

Now consider eq. (|C5|) in the limit that H > Z CR . Since we are primarily interested in the force near the surface, we 
will ignore the exp(—2H/l CR ) dependent terms since they only provide order unity corrections. We then find 

dp CR _ KQcr exp{-n) 
dz 3cZ CR 1 + y/e 

Dropping the y/e term, which is small, and using the above definition for Q CR , this simplifies to 

— — ~ — exp(-n). 

dz H A CR c 

We can see that force is reduce roughly exponentially with r„ as we move deeper in the atmosphere, eventually going 
to zero at the midplane. The maximum force is at the surface and is equal to 



^^~CR ^CR -^rnax 



dz H A CR c 



<X \ft. 



With the help of figures IC1I and IC2i the above expression has a clear interpretation. The force exerted by the CRs is 
still given by F CR /(X CR c), but F CR is reduced from the maximum flux -F max by a factor of l CR /H because only CRs 
produced within ~ l CR of the surface can contribute to F CR at the surface. As we move deeper within the atmosphere 
the force is reduced because downward moving CR from above increasingly offset the upward moving CRs from below. 
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APPENDIX D: MOMENTUM AND ENERGY CONSERVATION FOR OPTICALLY THICK 

RADIATION-DRIVEN WINDS 

Here, CRs are taken to provide the radiation pressure that ultimately drives out a given galaxy's gaseous component. 
By making the simplifying assumption that the force exerted by the CRs upon the gas greatly exceeds the gravitational 
and thermal pressure forces, the luminosity in CRs may be related to the outgoing momentum of the wind by 

P w = Mv w = -S3.T, (Dl) 

c 

where P w , M, and v w is the momentum, mass loss rate, and terminal velocity of the wind, respectively (Lamers & 
Cassinelli 1999). In eq. (|D1[) . the optical depth of the wind r cannot be arbitrarily large and is constrained by energy 
conservation of the flow i.e., 

^ch < \mv 2 w . (D2) 

It follows that Mv w < 2L CR /v w , which gives a upper limit for the optical depth 

t<— . (D3) 

In the extreme limit, where r = 2c/ v w , CRs are completely rcdshiftcd away due to adiabatic expansion since they 
have transferred all of their energy to the matter. 

APPENDIX E: TURBULENT POWER OF MOLECULAR CLOUDS REGULATED BY COSMIC RAYS 

The observation that the light to mass ratio 72 of both "fully populated" molecular clouds and bright dense regions of 
starbursting galaxies are equal imply the presence of a common feedback mechanism. Here, we explore the possibility 
of whether or not the production and transport of CRs within giant molecular clouds (GMCs) can regulate star 
formation. 

In the Galaxy, there is ~ 10 9 M Q of molecular gas responsible for the formation of stars and most of this mass is in 
the form of GMCs (Binney & Tremaine 1987). We work under the assumption that turbulence is the dominant form 
of pressure support and that the molecular cloud is in virial balance. 

The amount of power, or luminosity, required to maintain steady turbulence in a GMC, L tm b, may be expressed in 
terms of M c , vi c , and At tU rb; which represent the mass of the GMC, the turbulent velocity of the outer scale l c , and 
the eddy turnover time on the outer scale, respectively. We have 

L turb ~ Mc-r^~ Mc^ 2 x lO^Mewf ifo 1 cr § s~\ (El) 

where Ai tur b — l c / v i c and M$, w 3 and Z 10 is M c in units of 10 6 M Q , l c in units of 10 pc, vi c in units of 3kms _1 , 
respectively. 

Now, we estimate an upper limit for the amount of power in CRs, which may be injected into the GMC. As usual, 
the source of CR power is assumed to be core-collapse SNe and well as winds from massive stars. Again, we take into 
account that only ~ 10 — 20% of massive stars are embedded i.e., in physical contact with their natal GMCs, at any 
given moment (Wood & Churchwell 1989). From this, we infer that a similar fraction of the SNe and stellar wind 
power comes into contact with its ancestral GMC. Putting this together, the CR energy injection rate that a GMC is 
subject reads 

M 

£cr d - ?T7^ L c R * 6 x 10 36 eo.iM 6 Ad i9 e 6 m SP ergs" 1 . (E2) 

Here £ and M mo \ represent the fraction of massive stars that reside within their parent GMC at any given time and 
the entire molecular mass content of the Galaxy and furthermore, M mo i t g and £q.i are M mo \ and £ in units of 10 9 A/q 
and 0.1, respectively. 

Understanding the parameter £ deserves some attention. One interpertation is that the fraction £ results from the 
possibility that every massive star eventually leaves its molecular cloud of origin before it expires into a supernova. 
The other interpretation is that a fraction £ of massive stars remain embedded for their entire lifetime. With respect 
to the former, CR feedback in GMCs plays a role only if stellar winds from massive stars effectively produce CR 
protons. For the latter, CR feedback may be important if either core-collapse SNe or stellar are responsible for their 
production. Note that £ is likely larger in the latter case. 

On energetic grounds, the above expressions indicate that CRs indeed may be a crucial agent of self-regulation for 
GMCs. 
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